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1. Introduction 



Carrier dynamics in a biassed superlattice (SL) has remained an active topic during 
the last two decades. Knowledge of the underlying physical processes is necessary 
for a better understanding of perpendicular transport of carriers in multi-quantum 
well systems, as well as for successful development of applications, such as microwave 
radiation generation [TJ [2j [3l 0] , quantum computing OE] and matter manipulation on 
the particle level [7j. 

Studies focussed on coherent carrier dynamics have revealed some dynamical 
features of great interest (see, for example, [8] and references therein) arising from the 
interplay between interminiband oscillations and tunnelling to the continuum. These 
studies contribute to understanding of the famous quantum-mechanical problem of 
tunnelling in presence of dissipation PUDIIII] that has been considered for superlattices 
in pj]. 

Whereas typically such studies involve a number of approximations, we adopted 
a computational approach with few limitations beyond those implicit in the model of 
a single electron in a biassed periodic potential. This enabled us to consider resonant 
dynamics at high and moderate bias in a finite superlattice, specifically Rabi oscillations 
(RO) and resonant tunnelling (RT), and to establish a relation between these two 
fundamental types of interminiband transport by simple means. A carrier dynamics 
solution relying on wavepacket time evolution also brings out certain features that are 
difficult to calculate otherwise, such as carrier behaviour in the vicinity of a resonance, 
and dependence of the period of Rabi oscillations on resonance index. 

The main purpose of this paper is to further analyze and explain the results obtained 
using the methods of our earlier work [13]. In particular, our approach helps to study 
near-resonance behaviour in more detail and understand the damping mechanism of 
RO and ways to reduce it. It also contributes to a link between RO and RT, an 
instance where quantum transport theory has been lacking to date. These results can be 
applied to any system possessing a Wannier-Stark ladder structure including photonic 
crystals [H] and optical lattices [T5] . 



2. Resonant interminiband dynamics 

This work focusses on longitudinal motion of a single electron in a zero-temperature 
biassed superlattice. We solve the time-dependent Schrodinger equation (TDSE) for 
\& (x, t) in a biassed periodic potential 

oo 

V( x ) — X/ Vsl(x — nd) + xF , where 



Vn 



x + a/2 x — a/2 



V SL (x) = tanh — —!■ tanh '— Q(x) 0(d - x) . (1) 

2 L a a J 

where Q(x) is the Heaviside function vanishing for x < 0, F = —eE is a uniform bias, 

E is the electric field, and Vsl{ x ) is the model superlattice potential in a unit cell of 



Resonant carrier dynamics 



3 



the periodic system, with width d and barrier thickness a. 

We consider coherent electron transport and thus omit electron scattering and 
relaxation processes; numerous studies of carrier coherency limits have been done in 
the past [TBI EZ] • These destructive processes are weak enough at low temperatures, to 
allow up to 14 Rabi oscillations to occur [18] . 

Several GaAs/Gai-^Al^As layered heterostructures labelled X = A, B, C, D are 
considered, whose parameters are shown in Tableland their band structure in Figure [TJ 
The layer profile function Vsl avoids discontinuities in the potential, which allows 
easier programming of the solver, and is more realistic than the often used square 
barrier approximation. Our numerical solutions of Equation [1] use discrete transparent 
boundary conditions [19] , which reduces the size of the space domain in which we must 
operate. For details, see reference |13j . 

The norm of the quasibound part of the wavepacket p and miniband occupancy 
p v are convenient properties by which to monitor interminiband dynamics; miniband 
occupancy is the wavepacket projection onto a tight-binding miniband v: 

Pl/ = J2\(^(x,t)\w u k (x))\ 2 . 

k 

Here W*(x) stands for a Wannier-Stark (WS) quasibound state corresponding to energy 
level centered on the well with index k, and belonging to miniband v (i/—l,2, . . . ); 
our initial wavepacket is centered on the well with index 0. The tight-binding Wannier 
functions from miniband v are denoted as w u (x). 

To refer to an interminiband resonance originating from an anticrossing of energy 
levels E* and E* +n in biassed sample X, we will use the symbol 7?." (X) and we will 
denote the resonant bias (the value of bias at which the peak of a resonance is observed) 
as F n , n being the resonance index (n = 1,2, . . . ). When of little importance, some 
indices may be omitted for brevity. For convenience the symbol G will stand for inverse 
bias 1/F (G n = 1/F n ). Unless specified otherwise, time is measured in units of the 
Bloch period Tg = 2iih/Fd. 

It should be noted that computing a complete complex-energy spectrum for a 
biassed multilevel system is a significant task beyond the aims of the paper, as can be 
seen from previous works devoted to this problem (for example, the formalism in [20| [2Tj 
utilizes existence of poles of the system scattering matrix). We instead consider a 
simplified framework proven suitable for the purpose of studying transient transport 
processes. 

2.1. Rabi oscillations model 

Typically, Rabi Oscillations (RO) are a result of interminiband transitions under 
external radiation when its frequency approaches the system's intrinsic value uj\2 = 
(E2 — Ei)/h [22]. This pumping provides excitation of carriers and their subsequent 
spontaneous emission of photons, which is widely applied in quantum cascade lasers [23J. 
A comprehensive overview of RO in the two-miniband approximation has been given 
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Figure 1. Band structure of the superlattice potentials in absence of bias (from left 
to right: samples A, B, C and D). 

Name Vq, meV d, nm(ML) a, nm(ML) a, nm 



Sample A 
Sample B 
Sample C 
Sample D 



212 
250 
212 
212 



13.0 (46) 
17.3 (61) 
13.0 (46) 
13.0 (46) 



3.1 (11) 
2.5 (9) 
2.3 (8) 
2.5 (9) 



0.4 
0.4 
0.4 
0.4 



Table 1. Geometric parameters of the model potentials used in simulations. 
Barrier height of 212 meV corresponds to £=0.18 and of 250 meV to ir=0.3 in the 
GaAs/Gai-xAl^As structure; ML stands for monolayer thickness, 1 ML = 0.283 nm. 



by Ferreira and Bastard [21] and further mathematical details can be found in 
references [25l [261 127] . In a biassed SL, interminiband transitions may also occur in 
the absence of external radiation [281 [29], [30] : when energy levels from coupled Wannier- 
Stark ladders (WSL) align in neighbouring wells, a carrier can easily tunnel between 
them. 

For an undriven SL (i.e. not exposed to an external oscillating electric field) a 
simplified yet fruitful analogy to an unbiassed system under external radiation can be 
made. For illustration, consider Bloch oscillations of an electron in a biassed SL, where 
the electron oscillates in a Bloch domain, with turning points of the motion defined by 
the edges of a miniband. One can think of this process not as a motion in space, but 
as oscillation in the local kinetic energy between the upper and lower band edges. For 
example, let an electron at t — have the mean position (xq) and the fixed initial energy 
E = (Ejj + El)/2, half way between upper and lower edges, and let a positive bias F be 
applied in the positive x-direction. An electron moving towards larger values of x will 
eventually collide with the lower band edge El, as the latter rises with slope F. The 
amplitude of motion in space is L/2 = (E - E L )/F = (E v - E L )/(2F) = (6E/2)/F, 
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where SE is the miniband width and L is the Bloch Oscillation (BO) domain width. 
The two domains, space and energy, are directly related through the bias value F. 

Assuming very large carrier lifetime 1/T compared to a typical oscillation period 
for moderate bias, and taking the WS states to be stationary and orthogonal, we can 
apply the acceleration theorem h(dk/dt) = F, combined with the dispersion relation, 
to describe wavepacket evolution over a Bloch oscillation cycle. Making use of the 
decomposition V(x) = Vsl + Fx we can write 

When the lowest minibands are tightly bound (which is common for a semiconductor 
SL), carrier dispersion in MBn can be approximated as E(k) = E n + (5E/2) cos(kd+n7c) 
and the mean position evolves as (x(t)) — (x(t )} = (L/2) sin {unit — to)); where u B is 
the Bloch angular frequency. 

Thus, the local kinetic energy is now time- dependent: 

E(t — t ) = (E v + E L )/2 - F((x(t)) - (x(t ))] 

= E -(FL/2) sin (w B {t-t ) 



= E -(5E/2) sm(u} B (t-t )) . (3) 

The wavepacket's evolution in coordinate space at fixed energy can effectively be 
replaced by its evolution in kinetic energy space. With minor variations, the same 
viewpoint can be adopted for intrawell oscillations in a similar fashion. Note that the 
intrawell oscillation frequency is a natural frequency of the system: u nm = (E m — E n )/h. 

Given the fact that interminiband transitions of an electron are a result of quantum 
interference of Bloch and intrawell oscillations [13J, one can extend this analogy to 
Rabi oscillations. Then the SL would be subject to a field of natural frequency ou nm , 
modulated by the Bloch frequency u B . In presence of the ui nm harmonic, this system 
is poised to undergo an interminiband transition, with a certain Rabi frequency u B . 
This is similar to the well-understood system of a two-level atom under an external 
monochromatic radiation resonance field f n . The Schrodinger equation for such a system 
transforms into 

~ ih ^m^ = (^ + v)^(x,t) 

with V = x f n cos(u B t), (4) 

where Hq refers to an undriven biassed SL in the tight-binding approximation, x is the 
dipole transition operator between two energy levels (minibands in our case), and /„ is 
the amplitude of the external electric field. 

The factor xf n , interpreted as a transition operator in our effective model, 
represents the strength of interminiband coupling and depends on many factors, 
including the applied bias and the potential shape. For a small xf n , perturbation theory 
can be applied and an expression for the population of the second level evolving from an 
initial wavepacket ^(x, 0) = W\(x) can then be obtained in the same fashion as for an 
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irradiated two-level atom. In our notation, for a two-level system driven off its ground 
state [22] one obtains: 

nt 



P2 



(t) 



(— V 

V p ' 



C{G) sin" 



r pmax 
1 R 



P 

£{G) 

r 

rpmax 



P 



with 



i + ((G-G n )/r 



1 -1 



£(G) 



where 



a;on/(^2 - E x ) 

( WT(*) 



x I ^ 2 m ( 



(5) 



with n being the index of the resonance considered. These equations are 
straightforwardly extended to an arbitrary pair of interacting minibands, subject to 
validity of the perturbation theory used. It is interesting that in order to analytically 
predict the P2/ p{t) curve in the entire near- resonance region we require only the values 
of the dipole matrix element x nm and of {p2/p) max , computed at the resonant bias field. 

This simplistic derivation for BO in a two-level tightly bound system applies to the 
wavepacket's 'centre of mass'. We neglected any change in dispersion of the wavepacket 
in the process of BO and thus it is only an approximation. Since we consider moderate 
to high fields in this paper, carrier decay is rapid enough to reduce the influence of this 
factor. As will be seen in subsequent sections, this simple model provides a surprisingly 
good fit to the simulation data. 

In the case of a strong field, the domain width in space of the wavepacket's 'center- 
of-mass' BO is typically smaller than a potential cell width. This implies that: (i) 
energy levels are sparse due to large splitting of WSL; hence near a resonance there is 
one preferred tunnelling path between MB1 and MB2; (ii) at a resonance, WS states 
are reasonably localized and have an exponentially vanishing tail, so we expect that 
xo n oc e~ n (which has been predicted by a two-level atom model [29] and also has been 
explicitly calculated for a driven multiband SL [20]). 

The above argument assumed the system to be in a steady state. However, in 
practice an initial non-equilibrium configuration undergoes a relaxation process to a 
state with lower potential energy. We studied the evolution of an initial wavefunction 
ty(x,t = 0) = C\W\{x) + C2W2(x), a linear combination of tight-binding Wannier states. 
(Obtaining correct quasibound WS functions is not important for our primary goal to 
study steady-state RO dynamics.) The steady-state tunnelling process showed little 
dependence on the linear combination: components of which belong to higher 

minibands tunnel out rapidly during the initial relaxation period, which only scales 
down the norm of the quasibound wavepacket being observed. The coefficients c\ and 
C2 should be set so that after an initial period of relaxation of a non-equilibrium state, 
the resulting RO are clear and of sufficient magnitude. The maximum magnitude of RO 
between two lowest minibands reaches nearly unity for the two extreme cases |ci| = 1 
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and | C2 1 = 1, with less carrier amplitude decay for the former choice. Hence, we adopted 
the initial wavefunction ty(x, t = 0) = Wi(x) for most of our results. 

2.2. Near-resonance behaviour 

In a Wannier-Stark ladder, energy level anticrossings show a dramatic increase in carrier 
tunnelling rate, which is ascribed to resonant Zener tunnelling. When studied in 
detail, carrier dynamics at some of these anticrossings exhibits a prominent oscillatory 
pattern (Rabi oscillations) which decays rapidly. For clarity we will call the latter 
Rabi resonances and the others tunnelling resonances. While both are anticrossings of 
complex energy levels (complex due to their finite width), it will allow us to differentiate 
the two types. At strong fields, above-barrier tunnelling resonances start to occur; we 
will include them as tunnelling resonances also. 

A study of isolated Rabi resonances in [133 revealed that a carrier undergoing RO 
produces damped oscillations in miniband occupancy, of the form: 

p n /p (t) =P n + A R e^ At cos ( ^ + (f) n ) , 

with p (t) = exp ( - 7* - a (1 - e~ w ) ) (6) 

in agreement with the general theory [21]. Here Tr, Ar and 7^ are the period, amplitude 
and decay rate of Rabi oscillations respectively; <p n is an initial phase determined by a 
particular form of 0) {e.g. <p\ ~ for 0) = Wi(x)); 7 is the decay rate of the 
entire wavepacket and P n is the asymptotic value of p n /p(t) in the limit t — > 00. The 
term —a (1 — e _7ne ') represents a dip in the decay rate due to two-exponential decay [32] 
and was found to be vanishingly small close to the resonant bias. Its presence is due 
to the initial relaxation of wavepacket components orthogonal to eigenfunctions of the 
current miniband, with a being proportional to the norm of the orthogonal part and 7 ne 
being the relaxation rate, typically larger by two orders of magnitude than 7. This term 
compensates for the change in states' orthogonality at higher fields. A typical example 
of time evolution of near-resonant interminiband occupancy dynamics, and its fit using 




Figure 2. Occupancy functions (solid lines) and their fit using Equation [5] (broken 
lines almost completely hidden behind them) for ty(x,0) = wi(x) near a resonance 
1Zl 2 {A) at G = 0.44058 nm/meV (0.8 HWHM away from the resonant bias); the 
upper curve corresponds to n = 1 and the lower to n = 2. 
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Figure 3. Behaviour of key dynamical parameters for 1Z\ 2 {A) (left column) and 
TZ\ 2 (R) (right column). Top panel shows RO amplitude vs. inverse bias (chain-dotted 
line shows a Lorentzian fit C(G), solid line shows fit using Equation [9|) , second from 
top panel: relative RO amplitude vs. relative RO period (filled circles correspond 
to G < G3, empty circles to G > G3; solid line shows fit using Equation [5] given 
Equation [7|) , third from top panel: asymptotic occupancy of the first miniband vs. 
inverse bias (solid line shows fit using Equation [5]) , fourth from top panel: relative 
RO period vs. asymptotic occupancy of the first miniband (filled circles correspond 
to G < G3, empty circles to G > G3; solid line shows fit using Equation [5]) ; bottom 
panel: carrier decay rate vs. inverse bias. 
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Figure 4. Damping rate of RO vs. inverse bias for 1Zi 2 (A), in units of 1/Tr. 



Equation [6] is shown in Figure [2j 

Since the parameters in Equation [H] provide a good description of carrier dynamics 
for moderate bias, and we will investigate their behaviour near a resonance. It has been 
found previously [13] that the period Tr of Rabi oscillations shows a root-Lorentzian 
peak around the resonant bias: 



T R (G) = T™* yJC(G) , (7) 

with the Lorentzian C defined in Equation [5] We will refer to the parameter T of 
this Lorentzian as the HWHM (half-width at half-maximum) of the resonance under 
consideration. It has also been shown that the peak period of Rabi oscillations changes 
exponentially with resonance index n\ 



rpmax rpmax I rp 

1 n — 1 l \ lr < 



max I rpmax ' " ' 



For a double quantum well system, perturbation theory predicts that the frequency 
of Rabi oscillations strictly at resonant bias is given by Tuor/2 = \ V(x) \ *&r), 
the tunnelling matrix element through the barrier separating the two wells [31] and 
is essentially the splitting between the energy levels. Thus the inverse of T R nax of 
Equation [7] is the minimum energy level splitting at their anticrossing. 

According to these earlier findings, the occurrence of resonant bias values 
also correlates well with complex energy spectrum anticrossings. As an example, 
the bias values for anticrossings between two lowest minibands in sample A, 
obtained from our simulation (Gi = (6.9±0.2) nm/meV, (j2 = (3.4±0.2) nm/meV and 
G3=(2.33±0.05) nm/meV), are reasonably close to those calculated by K. Hino 
et al. [33] (7.2 nm/meV, 3.6 nm/meV and 2.4 nm/meV, respectively). This agreement 
demonstrates that interminiband motion prevails at complex energy anticrossings. 

The parameters Ar and P n derived from simulation results also have clear extrema. 
The following equations fit the data quite well (see Figure [3]): 



a a max 

A R - A R 



1 - C^(G) 

IT At, 

P n =pr x sm 



(8) 



sin 

OA max ' 
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n T™" 

Figure 5. Logarithmic fit of T n versus resonance index for 7£i2(A) (left) and its 
relation to T% ax (right). 

with and P™ ax being the peak values; P n = P n for the lower resonantly coupled 

MB and 1 — P n for the upper one. The other two parameters, 7 and 7,4, are extremely 
sensitive to coupling to higher minibands and their bias detuning dependence varies 
from one resonance to another. However, at a resonant bias, 7,4 always reached its 
virtually zero minimum, and 7 its maximum (see Figures [U 0J) . 




0.235 0.24 0.245 0.25 0.235 0.24 0.245 

G (nm/meV) G (nm/meV) 

Figure 6. Decay rate 7 of a wavepacket vs. inverse bias, near TZf 3 (C) (left) and 
fc? 3 (D) (right). ' 




G (nm/meV) G (nm/meV) G (nm/meV) 



Figure 7. Near- resonant dynamics of second miniband occupancy at 7?.f 2 (A) (left), 
n\ 2 {C) (center) and 7^ 2 (B) (right). 
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G-G 2 (x1 nm/neV) G-G 2 (x1 nm/|ieV) 

Figure 8. Amplitude of RO and of first miniband occupancy near resonances 1Z\ 2 
and 1Z\ 3 from middle panel of Figure [7] Solid lines show a superposition of Lorentzian 
curves best fitting simulation data, in circles; broken line is an estimate of a fit in the 
absence of 7£?o. 



As the index n rises, the TZ^ are expected to narrow, since the length of the 
tunnelling pathway into the adjacent miniband increases and so does the system's 
sensitivity to bias detuning, due to related frequency detuning of the Bloch oscillations; 
this is reflected in a decrease of x 0n in Equation [51 As predicted by Equation [SJ the 
resonance HWHM dependence on n closely follows an exponential law (Figure [5]) with 

r n = r\ (iyr^ 1 (9) 

The anticipated relation 

rr x = G /T n (10) 

also holds quite well: fitting Equation (Q produced a value T™ aa T n =(139 ± 5) nm/ /ieV 
compared to Go =(145 ± 1) nm/ /zeV for sample A, and calculations for sample B showed 
that T 3 maa T 3 =(1.05 ± 0.01) nm/meV compared to G =(1.02 ± 0.01) nm/meV. 

When an ensemble of resonantly coupled minibands is weakly bound, Rabi 
oscillations are weak and overdamped and it is mostly resonant tunnelling that is 
seen [M| |35| [36] to prevail over RO; these represent tunnelling resonances. In this case 
the carrier wavepacket escapes to the continuum very quickly and exhibits no persistent 
RO. Then 7 becomes the key parameter in the dynamical description and features a 
clear peak centered at the resonant bias (Figure EJ). Sensitivity of 7 to coupling to higher 
minibands makes the shape of curves 7(G) vary for different resonances. 

While we have used perturbation theory to interpret our calculated results, the good 
fits obtained by using Equations ( l9lfT0]) suggest that the model remains valid at strong 
biases as well. We attribute this to the fact that the structure of the Wannier-Stark 
ladder is preserved over the bias domain considered [13] ■ At stronger electric fields, 
however, Bloch oscillations are replaced by sequential tunnelling, and projection on 
minibands refers merely to the wavepacket's distribution in energy (and hence between 
wells in real space), rather than a decomposition into Wannier-Stark states. Indeed, the 
n th tight-bound miniband contains only harmonics with wavelengths A G [|g?, ^-d]. 
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3. Superposition of carrier dynamical patterns 

So far we have considered isolated interminiband resonances. They occur typically 
in a strong potential at moderate fields where resonances have narrow HWHM 
compared to their spacing in inverse bias space, and hence they rarely overlap. It 
is instructive, however, to consider higher-field resonances in a weaker potential (e.g. 
Alo.3Gao.7As, sample B in Table 1), where the interaction between the two resonantly 
coupled minibands represents only one of several prominent interference paths. Under 
these conditions, the interaction between these paths becomes important and largely 
determines the overall dynamics. This gives rise to some particular dynamical patterns 
that we interpret as interference of tunnelling and Rabi resonances, between different 
pairs of minibands situated in close proximity. 

3.1. Closely situated resonances 

We now consider a case of strongly interfering paths, to illustrate how different dynamics 
contribute in such a mixed mode, and gain understanding of RO stability. This situation 
typically occurs when the HWHM of one resonance is large compared to the HWHM of 
another resonance that is coupled to the higher minibands. Then one dynamical pattern 
is seen to be superimposed on another, as in the centre panel of Figure [7] where one 
can see evidence of both oscillatory and tunnelling carrier behaviour. The oscillatory 
pattern related to the 1Z\ 2 resonance, vanishes at the 7£f 3 peak, and the Pi(G) and 
Ar(G) curves in Figure [8] demonstrate a sharp extremum resembling a superposition of 
two Lorentzian-like curves. 

It is worth mentioning that, in the centre panel of Figured the tunnelling resonance 
7£f 3 that reaches across three minibands, lies above the barrier height of the SL. 
This gives it comparable strength to 1Z\ 2 acting between adjacent minibands, which 
significantly alters the Rabi resonance pattern. Clearly, resonances which cross three 
minibands are much weaker than those that cross just two, and a necessary condition 
for the former to affect the latter significantly, is strong bias resulting in weak carrier 




G 



Figure 9. Schematic drawing to illustrate the mutual impact of energy level 
anticrossings. Chain-dotted line shows position of level E® in the absence of repulsion 
from level £3, shown as a thin solid line. For further explanation, see text. 
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localization. Despite such a strong field, we can use a Wannier-Stark ladder to describe 
transitions of an electron for the strong fields considered in Figure [7J since the level 
structure disrupts only at fields high enough for the potential drop per unit cell to be 
close to the smallest interminiband separation in the system E 2 — E\ [13]. In our case, 
this value corresponds to the threshold inverse field G = 0.28 nm/meV in sample B (the 
strongest potential considered with smallest E 2 — Ei separation). 

In Figure [7J we present a few typical cases of interminiband resonance occurrences. 
The leftmost panel is an example of an isolated Rabi resonance having an almost 
perfectly symmetric structure; the two panels to the right are examples of interference 
between Rabi and tunnelling resonances. In the centre panel, the tunnelling resonance 
1Z\ Z appears to be particularly strong and dominating. Hence, to the right of the Rabi 
resonance 1Z\ 2 peak, one sees a clear oscillatory pattern, whereas to the left of the peak 
it is smeared out by strong tunnelling. The close presence of the tunnelling resonance 
1Z\ 3 has also the effect of shifting the Rabi resonance 1Z\ 2 peak along the bias scale. 

From the relation F n = uFq, the resonance 1Z\ 2 (centre panel of Figure CJ) should 
have a peak at G — (0.284 ± 0.002) nm/meV; however as a fit to the maximum period 
of oscillations indicates, it occurs at a lower bias G = (0.295 ± 0.002) nm/meV. The 
shift occurs because at the energy level anticrossing corresponding to 7£f 3 , the levels Eq 
and fif are repelled from each other [20J. At the same time, the position of level E\ 
involved in a Rabi resonance remains almost unchanged. Provided that the minimum 
mismatch of coupled energy levels {E\ — Eq) produces the longest period of RO [13J, 
this results in a shift of the Rabi resonance 1Z\ 2 peak as qualitatively shown in Figure [9j 
There the minimum energy difference at 1Z\ 2 , AE [ d , located at G Q u in the absence of 
a lZl 3 anticrossing, shifts to G new due to superposition of the energy level anticrossings 
K\ 2 and K\ 3 . Since the repulsion from the two anticrossings adds constructively, AE new 
is larger than AE id, and since the RO period is inversely proportional to the mismatch 
in energy level alignment, the peak period of RO shifts to a new value of bias. Thus 
the proximity of another resonance has a twofold effect on the resonance considered: 
it shifts its resonant bias and reduces the maximum period of RO. The latter effect 
is demonstrated in the right panel of Figure [7J where a weak and narrow tunnelling 
resonance spread over G = 0.322 . . . 0.330 nm/meV reduces the RO period. 

3.2. Conditions for Rabi resonance 

The simulation data reveal that the decay of RO is slower far from tunnelling resonances 
and in stronger potentials, with tightly bound lowest minibands reducing the carrier 
decay rate 7. At the same time, the RO decay rate nearly vanishes at the very peak of a 
Rabi resonance, where the wavepacket decay rate is at its highest. Thus tunnelling rate 
alone is not a reliable indicator of Rabi oscillation damping. As Figure fTUl demonstrates, 
resonant tunnelling and Rabi oscillation patterns may occur for the same wavepacket 
decay rate. 

In the centre panel of Figure [7] at inverse bias G^ = 0.312 nm/meV there is a 
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clear RO pattern, whereas at the bias G^ = 0.274 nm/meV (which is symmetric to 
j n relation to the Rabi resonance peak located at G — 0.293 nm/meV), Rabi 
oscillations die out quickly. This behaviour goes beyond the symmetric structure of 
an isolated resonance as in Equation [7J (left panel of Figure [7]); it is interference of a 
tunnelling resonance at G^ that overdamps the Rabi oscillations. To generalize, the 
alignment of energy levels from higher minibands sets the wavepacket behaviour model 
and the link between RO and tunnelling resonance (both being a product of wavepacket 
self-interference) is determined by a particular arrangement of Wannier-Stark ladders: 
specifically, by the ratio of interminiband tunnelling rates, as will be argued next. 

For illustrative purposes, we turn again to the two-miniband model and denote 
the rate of tunnelling from MBz/ to MB// as 7„ M ; this is a convenient measure of 
strength of interminiband coupling. The magnitude of 7^ is inversely proportional 
to the difference in energy between the resonantly coupled energy levels from WSLz/ 
and WSL/z [37J. In the centre panel of Figure [7J these level pairs are: E® and E\ for 
the Rabi resonance and E\ and E\ for the tunnelling resonance. As seen in Figure [TUl 
for the inverse biases G^ and G^ featuring the same carrier decay rate, the difference 
(E% — E®) — AE21 determining the tunnelling between MB1 and MB2 is -2 meV and 
6.2 meV, respectively (AE 2 i is the energy spacing between the centres of MB2 and 
MB1, which equals 88.8 meV for potential A), whereas at the same values of bias, 
{El - El) - AE 31 for the ensemble of MB1 and MB3 equals -10 meV and 10 meV. 
Considering that the widths of the tight-binding MB1, MB2 and MB3 are 8.8, 39.4 and 
98.6 meV, the disparity between the two values of [E\ — E®) — AE 2 \ oc I/712 is much 
more significant than that of (E^ — E®) — AE 3 i oc I/713. Hence, the ratio 712/723, which 
indicates the isolation of the coupled bands MB1 and MB2 from higher minibands, is 
much smaller in the case of G^\ Due to a weaker isolation of the ensemble at G^\ 
the quasibound oscillating part of the wavepacket is more subject to tunnelling into the 
continuum, so RO are destroyed and turn more quickly into resonant tunnelling. The 
evidence for this is the middle line in Figure [10] featuring strongly damped oscillations. 
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Figure 10. Dynamics of first miniband occupancy in the vicinity of lZf 2 , from the 
centre panel of Figure[7]at G^ — 0.301 nm/meV (upper curve), G^ = 0.293 nm/meV 
(middle curve) and G^ = 0.274 nm/meV (lower curve). The upper [lower] curve is 
shifted upwards [downwards] by 0.1 for visibility. 
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Figure 11. Ratio of interminiband tunnelling rates 712/723 vs. inverse bias (in 
resonance HWHM units), from Equation ITTI given Equation [9l Solid and broken lines 
show the cases when Pi = 0.5 and 0.6 at resonance (ideal and nearly ideal cases), 
respectively. 



This also explains why, at the peak of an isolated Rabi resonance, the RO damping 
rate is the lowest. For example, at R12 two lowest minibands merge through resonant 
coupling, so 7i 2 and 7 are exploding while 7 23 changes relatively little. In other words, 
at an anticrossing, the decay rate of a more stable state changes much more than that of 
a less stable one. Thus, a surging ratio of 712/723 makes RO decay go nearly to zero at 
the peak, despite homogeneous level broadening which reaches its greatest extent there. 
Away from the peak 7 falls off and RO persistence decreases, because of the quick drop 
in 712 (Figure H]). This is a somewhat counterintuitive example of how a quasibound 
state with shorter lifetime can produce more prominent oscillations. 

Since 712 is closely related to the interminiband transition matrix element 
V12, we can estimate the near-resonant behaviour of the latter by invoking some 
dynamical equilibrium considerations. From an elementary two-miniband model (see 
the Appendix) we obtain 

712 2Pi - 1 

— = (11) 

723 Pi(l - Pi) 1 ] 

Neglecting tunnelling pathways between poorly aligned energy levels (e.g. E\ -> E\ 

at 7^12) and over many potential wells (in case of resonance across three minibands), 

V12 = 712 ■ Assuming that 723 changes over bias near TZ\2 adiabatically slowly compared 

to 712, we get V\% oc p^ff-p^ ; with the fit for P\{G) from Equation^ we can estimate 

the behaviour of Vu as shown in Figure [TT1 In the ideal case of Pi = 0.5 at the peak, 712 

vanishes, which reflects a complete merger of the two lowest minibands. This example 

shows a link between the results obtained for bias detuning dependence of the key 

dynamical parameters, and theoretical values. With few simplifying approximations, 

one can investigate near-resonant behaviour of theoretical parameters based on carrier 

dynamics in a finite superlattice system in the above manner, using a system of coupled 

damped oscillators described in [38] . 
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3.3. RO damping and quantum interference 

On a fundamental level, RO damping may be seen to originate from self-interference of 
a wavepacket that is facilitated by interference between small fractions of wavefunction 
that a wavepacket constantly emits through intrawell and Bloch oscillations (we will 
call them 'ripples'). A wavepacket initially set in the ground miniband of a biassed SL, 
starts to leak out into the second miniband, since the interference of the 'ripples' is 
constructive at the location of the initially unpopulated MB2, and destructive at MB1. 
When most of the wavepacket has tunneled out into MB2, the tunnelling direction 
reverses. This follows since by then most of the emission is from MB2. In the ideal case 
of zero RO damping, this reversion is a mirror process; in practice, the tunnelling rate 
to higher minibands cannot be neglected. This causes the amount of emitted 'ripples' 
to reduce faster in MB2; such an inequality breaks the mirror symmetry of the two 
transitions. The oscillation pattern smears out and RO are damped more quickly, the 
larger is the imbalance in the rate of 'ripples' escaping to the continuum. That is, 
inversely proportional to 712/723- In other words, if 723 cannot be neglected compared 
to 7 12 , it breaks the anisotropy in the system's tunnelling pathways: in addition to 
MB1^MB2, there appear two additional pathways, MB1^MB3 and MB2<->MB3 with 
very different tunnelling rates. 

4. Conclusion 

From the results of numerical simulations of near-resonant carrier dynamics at an 
isolated resonance, we have proposed a set of equations governing wavepacket behaviour; 
the dynamical parameters near a resonance were found to exhibit extrema of various 
shapes. For overlapping resonances, a superposition of energy level anticrossings 
produces a shift in resonant bias. It also reduces the maximum period of RO and 
perturbs the dependence of dynamical parameters on bias. A superposition of Rabi 
and tunnelling resonances also allowed us to observe the transition between oscillatory 
and tunnelling coherent wavepacket dynamics, and to examine the mechanism of RO 
damping. 

Persistence of Rabi oscillations near a resonance was demonstrated to be 
independent of homogeneous level broadening and to depend on the ratio of 
interminiband tunnelling rates 712/7235 which can serve as an estimate of the RO 
damping rate. 
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Appendix: dynamical equilibrium for a two-miniband system 

We consider an elementary system of two strongly interacting minibands coupled to 
continuum states. Further, let us denote the tunnelling rate from MBz/ to MB/i as 7^, 
from MBz/ to continuum as 7^ and the current relative occupancy of MBz/ as p v j p = p v . 
In steady tunnelling mode, the miniband occupancies are in dynamical equilibrium. 
Neglecting direct tunnelling from MB1 to the continuum, the flow of probability is 
balanced as follows: 



where p has components pi, p 2 , and the second equation expresses the condition for 
dynamical equilibrium. A little algebra produces the quadratic equation 



with k — p~\j pi > 1. Taking p 2 = 1 — p\ and assuming j 12 ~ 721 = 7 (which holds well 
in the vicinity of a resonance peak), in the steady-state limit i^oowe arrive at 



with the notation a = 700/27 and 700 = 7200- Note that the extreme case of RT 
corresponds to a — > 00 with P 1 — > 1, and of RO to a — > with P 1 — > 0.5. 
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